Conductance of ion channels and nanopores with charged walls: a toy model 
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We consider ion transport through protein ion channels in lipid membranes and water-filled 
nanopores in silicon films. It is known that, due to the large ratio of dielectric constants of water 
and the surrounding material, an ion placed inside the channel faces a large electrostatic self-energy 
barrier. The barrier leads to an exponentially large resistance of the channel. We study reduction 
of the electrostatic barrier by immobile charges located on the internal walls of the channel. We 
show that the barrier practically vanishes already at relatively small concentration of wall charges. 



Protein ion channels functioning in biological lipid 
membranes is a major frontier of biophysics An ion 
channel can be inserted in an artificial membrane in vitro 
and studied by physical methods. Similar artificial de- 
vices - water-filled nanopores, are studied in silicon, sili- 
con oxide films and polymer membranes^. In both cases, 
one can study a single water filled channel connecting 
two reservoirs with salty water (Fig.^l. A static voltage 
applied between these reservoirs drops almost entirely in 
the channel due to the high conductivity of the bulk solu- 
tion. The voltage drives salt cations and anions through 
the channel. One can measure the ohmic resistance of 
the channel. 

This resistance may be exponentially large due to the 
fact that the dielectric constant of water m ~ 80 greatly 
exceeds that of the surrounding media («2 — 2 for 
lipids and k,i ~ 4 for silicon oxide). Indeed, in this case 
the electric field of an ion traversing the channel is forced 
to stay inside the channel (Fig.GJ. This creates a barrier 
U(x), where a; € [— L/2, L/2] is the ion coordinate inside 
the channel. The barrier is the difference between the 
self-energy of the ion at the point x inside the channel 
and the self- energy in the bulk 4 . It is the maximum of 
the barrier, [7(0) = Ul, that determines the resistance of 
the channel. If a channel is very long the electric lines 
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FIG. 1: Electric field of a cation in a short cylindrical channel 
with the large dielectric constant ki ^> «2- L is the channel 
length, a is its radius. The self-energy barrier is shown as a 
function of the coordinate x. 

leak through the protein walls and lipids so that the bar- 
rier saturates as a function of L. Roughly speaking this 
happens 4 *^ at L ~ a(ni/ , where a is the chan- 
nel radius. In this paper we assume that the channel is 



shorter, so that we can neglect the field leakage. In this 
case, calculation of the barrier height is very simple^. 
The electric field at a distance x > a from a cation lo- 
cated in the middle of the short channel is uniform and 
and according to the Gauss theorem is Eq — 2e/nia 2 . 
The energy of such field in the volume of the channel is 
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where the zero argument is added to indicate that there 
are no other charges in the channel. Ul (0) is proportional 
to L and (for a narrow channel) can be much larger than 
ksT, making the channel resistance exponentially large. 

At large concentration of salt in surrounding water the 
electrostatic barrier is reduced by screening^!. In bio- 
logical channels the nature uses more effective approach. 
Channels designed for the transport of cations (K, Na, 
Ca) have negative charges on internal walls. For exam- 
ple, the potassium channel has 8 amino-acids with neg- 
atively charged radicals build into the wall of the pro- 
teinic. Walls of artificial nanopores, generally speaking, 
are charged as well and one can control these charges by 
a chemical treatment and/or tuning pH of the solution. 
The goal of this paper is to study the effect of immobile 
wall charges on the electrostatic barrier and the channel 
resistance. For certainty we assume that wall charges are 
negative and equally spaced along the channel with the 
linear density n w . 

We show below that in a large range of salt concentra- 
tions, c, the wall charges attract equal number of cations 
from the solution in order to make the channel neutral. 
Our theory is based on the observation that the Coulomb 
interaction of all charges in the short channel obeys the 
one dimensional Coulomb law: <fr(x) ~ \x\, same as for 
parallel uniformly charged planes. Indeed, let us consider 
a negative charge fixed at the wall and a cation, which 
arrived to the channel in order to screen it (Fig. [2} . The 
uniform electric field between them creates the confin- 
ing "string" potential Q(x) — eEo\x\. This situation 
reminds two quarks confined in a meson. Condition 
Q(x) = ksT defines the characteristic thermal length 
of such classical "atom", xt — ksT/eEo — a 2 /Is, where 
Ib = e 2 / (kiUbT) is the Bjerrum length (for water at the 
room temperature Ib — 0.7 nm). This "atom" is similar 
to an acceptor in a semiconductor (the classical length 
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xt plays the role of the effective acceptor Bohr radius). 
It is clear that at a small dimensionless concentration of 
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FIG. 2: A cation bound to a negative wall charge (circled). 
When the cation moves away from the host the energy grows 
linearly with the separation x. 

wall charges 7 = u w xt <C 1, each of them binds only 
one cation (Fig. [3^). Resulting neutral atoms do not in- 
teract or overlap with each other. This system reminds 
a lightly doped p-type semiconductor at very low tem- 
peratures when all holes are located at their acceptors. 
Let us show that already for relatively small 7 < 1 the 
electrostatic barrier U l may be substantially reduced. 




state. This difference is reduced by a large entropy S of 
the saddle point configurations. 

To simplify the calculations let us imagine that the wall 
charges form a periodic one-dimensional lattice along the 
x-axis. Then F s = Ul{0) — TS, where the entropy of 
the channel, enhanced by the charge unbinding, is S = 
n w L ln(ira 2 /2n w v). This gives for 7 <C 1: 

Ul{i) = F S -F g = U L (0)[1 - 4 7 ln(l/2 7 )] . (2) 

In the opposite limit, 7 > 1, one may expect that 
atoms overlap and destroy each other making cations 
free. In other words, one could expect an insulator- 
to-metal (or deconfinement) transition at some critical 
7 ~ 1. This does not happen, however. We show be- 
low that, due to the peculiar nature of the Id Coulomb 
potential, the barrier proportional to the system's length 
persists to any concentration of the wall charges, no mat- 
ter how large it is. Its magnitude, though, decreases ex- 
ponentially at 7 3> 1, 

U L (rf)/U L (p) cx exp(-11.03V7) ■ (3) 

Using numerical procedure outlined in the end of the pa- 
per we calculated the ratio 7(7) = Ul(j)/Ul(0) at any 
7 and plotted it along with the asymptotic Eq. J5J in 
Fig. 0] Let us discuss the range of the salt concentration 



FIG. 3: The ground state and the transport saddle point of 
the channel. Only right half of the channel is shown, (a) The 
ground state: all cations are bound to wall charges (circled), 
(b) The transport saddle point: an extra cation in the middle 
of the channel (left) makes other cations free. 

The equilibrium partition function of the channel is 
Z = (2ira 2 x T /v) n ™ L , where v is the normalization cell 
volume. This gives for free energy in the ground state 
F g = —4UL(0)j\n(2ira 2 XT/v). To evaluate the transport 
barrier one needs to know the free energy conditioned to 
the situation where an extra cation is placed in the mid- 
dle of the channel, F s . It creates the electric field Eq 
which orients the dipole moments of all "atoms" along 
its direction. In other words, it orders all charges in an 
alternating sequence of positive and negative ones. Due 
to the Id nature of the problem, this field unbinds each 
cation from its wall host and makes it free to move be- 
tween nearest neighbor wall charges (Fig. Eb). Indeed, 
according to the Gauss theorem the wall charge closest 
to the extra cation changes electric field from Eq to —Eq, 
then its cation changes it back to Eq and so on. Thus, 
at any position of cations between their nearest neighbor 
wall charges, the electric field is E = ±Eq, i. e. \E\ 
is constant throughout the channel. Therefore, the to- 
tal electrostatic energy is again given by Eq. (JIJ, One 
could think that the transport barrier is still given by 
Ul(0). This is incorrect, because the barrier is actually 
determined by the difference of the free energies F s — F g 
of the collective transport saddle point and the ground 




FIG. 4: The function /( 7 ) = U L (j)/U L (0). Its 7 < 1 asymp- 
totics, Eq. JSJ, is shown by the dotted line. The inset shows 
how asymptotics of Eq. ||3J (dashed line) is approached at 
7> 1. 

in the bulk solution c, where the above results are valid. 
A convenient dimensionless variable for c is a = irca 2 XT- 
Eq. J2| is valid only when 7>a, when due to neutrality 
the total number of cations in the channel is close to n w L. 
In the opposite case, a > 7, additional cations together 
with equal number of anions enter the channel. In order 
to minimize the energy all positive and negative charges 
should alternate along the length of the channel or in 
other words they have to be ordered. This means that 
each segment of the channel between two nearest neigh- 
bor wall charges gets an integer number k of additional 
cation-anion pairs. 

Let us now calculate the linear in L transport barrier 
t/i(a, 7) at arbitrary a, 7 <C 1. At the transport saddle 
point, when an extra cation resides in the center of chan- 
nel and creates the field Eq in both directions, all mobile 
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charges become free to move while keeping alternating in 
charge order. The energy of all such states is again equal 
to Ul(0), because l^l = E everywhere. This barrier is 
reduced by a significant entropy term. To evaluate it let 
us first consider entropy of a segment between two wall 
charges. If k cation-anions pairs reside in the segment, 
the total number of free ions there is 2k + 1. The total 
entropy of these 2k + 1 ions is 



k B ln[(7ra 2 /n w v) 2k+1 /{2k + 1)!]. 



(4) 



The factor l/(2k+ 1)1 reflects the fact that when cations 
and anions are ordered it is impossible to obtain a new 
configuration by exchanging them. Moving 2k ions from 
the bulk leads to the entropy loss 2kk B ln(l/cu) there. 
Subtracting this entropy from Eq. and taking sum 
over all k from to oo we obtain the total entropy of the 
saddle point configuration and the renormalized barrier 
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In the case a <C 7 Eq. JSJ matches Eq. (J2J. In the 
opposite case a » 7 Eq. JSJ crosses over to the result 
Ul{o) = Ul(0)(1 — 4a) obtained previously for an un- 
charged channel 6 -. 

So far we presented results for a periodic lattice of wall 
charges. To understand the role of random distribution 
of wall charges along the x-axis let us return to the case 
a < 7 C 1. It is easy to show that averaging over ran- 
dom nearest neighbor distances substitutes Eq. (J2J by 
U L {i) = U L (0)[1 - 4 7 ln(e- c /27)], where C = 0.577 is 
the Euler constant. Thus, the result for randomly dis- 
tributed wall charges is similar to those for the periodic 
one. 

Until now we concentrated on the barrier proportional 
to the channel length L. If a <C 7 there is an additional, 
independent on L, contribution to the transport barrier. 
It is related to a large difference of concentrations of 
cation inside and outside the channel. Corresponding 
contact (Donnan) potential Ud is created by double lay- 
ers at each end consisting of one or more negative wall 
charges and screening (positive) charge in water. 

For 7 <C 1 one finds \Ud\ <C Ul(j) and the channel 
resistance remains exponentially large. When 7 grows 
the barrier Ul{i) decreases and becomes smaller than 
Ud = ~ ferin(7/a), which increases with 7 and makes 
the channel strongly cation selective. In this case the 
measured resistance may be even smaller than the naive 
geometrical diffusion resistance of the channel. 

Let us, for example, consider a channel with L = 5 
nm, a — 0.7 nm, xt — 0.35 nm at c = 0.1 M and n w = 1 
nm -1 (5 wall charges in the channel), which corresponds 
to a = 0.035 and 7 = 0.35. The bare barrier Ul{0) — 
3.5k B T is reduced down to U L (^) = 0.2k B T. At the 
same time Ud = — 2.bk B T. Thus due to 5 wall charges, 
instead of the bare parabolic barrier of Fig. ^ we arrived 
at the wide well with the almost flat bottom (Fig.[5J). 

The contact potential Ud may be augmented by the 
negative surface charge of the lipid membrane^ or by 
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FIG. 5: The electrostatic potential for cations for the channel 
with 5 wall charges considered in the text. 



affinity of internal walls to a selected ion, due to ion- 
specific short range interactions^. It seems that biolog- 
ical channels have evolved to compensate large electro- 
static barrier by combined effect of Ud and short range 
potentials. Our theory is helpful if one wants to study 
different components of the barrier or modify a channel. 
In narrow artificial nanopores there is no reason for com- 
pensation of electrostatic barrier. In this case, our theory 
may be verified by titration of wall charges. Nanopores 
can also be longer than cl{k\/ kz) 1 / 2 so that electric field 
lines leakage through the walls becomes substantial. This 
leads to flattering of the parabolic barrier, but its depen- 
dencies on a and 7 remain qualitatively the same 6 . 

Let us elaborate now on the technical aspects of the 
derivations. As was first realized in Ref. Q, the parti- 
tion function of the Id mobile Coulomb plasma may be 
written as a trace of an (imaginary time) "evolution" op- 
erator, Z(q) = Tr cxp{— H (q)L/xT} with the Mathieu 
Hamiltonian H(q) = (ide — q) 2 — 2a cos 9. The variable 
q has a meaning of a not-necessarily-integer screening 
charge induced at the channel opening 6 . The ground 
state of the channel corresponds to q = 0, while the col- 
lective saddle point, having an uncompensated cation in 
the middle of the channel, to q = 1/2. The transport 
barrier is given by the difference between the two free 
energies: F s (l/2) — F g (0). As a result, the transport 
barrier of an uncharged channel is proportional to the 
width of the lowest Mathieu-Bloch band. It is a rapidl 
decreasing function of the mobile salt concentration a 

The immobile charges are represented by the charge 
"creation" operators exp{±i#}, where the sign is given by 
the sign of the static unit charge. For example, the par- 
tition function of the channel with the negative charges 
fixed at the positions x±,X2,--- is given by Z(q) = 

Tr ^ e -H(x 1 +L/2)/x T e -i6 e -H{x 2 -x 1 )/x T e -i0 _ _ \ j n 

simplest case of periodically placed charges with the di- 
mensionless concentration 7 one faces the spectral prob- 
lem for the non-hermitian operator U (q) — er* 1 ^! 1 e~ % ® . 
The barrier may be expressed in terms of its largest eigen- 
value e~ A 9 ' as Ul(j) = k B T(\fJ 2 — X^)jL/xt, where 

U{q)^ n = e~ A 9 '^n. One may demonstrate that, despite 
of being non-hermitian, the operator IA (q) possesses only 
real eigenvalues. 

In the limit of small concentration, 7 < 1, one may 
write the operator hi in the eigenbasis of the Hamiltonian 
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H (Bloch basis). Due to 7 < 1 condition only lowest 
eigenvalues of H should be retained. At q = (the band 
minimum) there is a well-separated ground-state of H, 
given by e °^ ~ 0. Evaluating the matrix element of 

e~ ld , one obtains Aq '' = — ln(2a). On the other hand, 
at q = 1/2 (the band maximum) and a C 1 there are 
two almost degenerate eigenvalues at cfi'2 — 1/4 =F a. 
Therefore, Xy 2 — 1/(47) — ln(sinha/7). As a result, 
one arrives at Eq. which in the limit a <C 7 yields 
Eq. ©. 

In the large concentration limit, 7 > 1, one may em- 
ploy a variant of the WKB approximation to find a spec- 
trum of Li. To this end one writes the eigenfunction as 
^(9) = e l ^ s ( s > and retains only the leading order in 7. 
The corresponding equation for S(9) reads as: 

[S'(d)] 2 +id-X=-(e ie + e- te ) . (6) 
7 

The two terms on its right hand side represent positive 
and negative mobile ions correspondingly. For large con- 
centration of fixed negative charges, 7 > 1, and moderate 
concentration of mobile salt, a < 7, one may consider 
only positive ions entering the channel and disregard the 
negative ones. This amounts omitting the e~ lS term on 
the r.h.s. of Eq. © (the fixed negative charge is encoded 
in the i9 term on the l.h.s.). Then shifting the 9 variable 
in the complex plane: 9^9 + Unfa /a) and defining 
A = A — ln(7/a), one brings Eq. Jfjjl into the parameter- 
less form: S'(9) = vA — id + e lB . The right hand side 
of this expression is an analog of the canonical momenta 
in the hermitian WKB scheme. In particular its zeros 
play the role of the classical turning points and deter- 
mine the structure of the branch cuts on the complex 
6*-plane. The integrals of the canonical momenta along 
such brunch cuts determine both the spectrum and the 
band-width upon changing the boundary parameter q. 

In the strip of the complex 6>-plane bounded by |9?0| < 
7r there are two turning points at 9 = ±6 — i<5 2 /6, where 



S <C 1 is defined as A = — 1 + 5 2 /2. The structure is then 
periodically replicated outside this strip. Choosing the 
brunch cut to run between the two turning points and ap- 
plying the Bohr-Sommcrfcld quantization rule, one finds 
the spectrum: A (n) = -l + (n + l/2) ^2/7. Since A -> -1 
in the limit 7 > 1, the two turning points are essen- 
tially close to the origin. To find the band-width one 
needs the "tunnelling" probability between the adjacent 
strips of the 6*-plane. To this end one chooses the branch 
cuts, which emanates from the turning points outwards, 
approaching 9 = ±n — ioo. The (exponentiated) inte- 
gral of the canonical momenta along such a cut gives the 
WKB tunnelling probability and hence the band-width 
Xf'L — X^K Such an integration leads to Eq. ©. 

For numerical calculation it is convenient to choose 
the basis of the angular momentum, e lme , to evalu- 
ate Z(q). In this basis the charge creation operator 
e~ 10 takes the matrix form [<5 mjm '+i], while the Hamilto- 
nian H„ hm < = [(to + q) 2 S m , m > - aS m , m '+i - aSm,m'-i\ ■ 
Truncating these infinite matrices with some large cut- 
off, one may directly exponentiate, multiply and trace 
them to find the free energy for any arrangement of fixed 
charges. 

As emphasized by the title of this paper we study here 
only a very simple model of a channel with charged walls. 
This is the price for asymptotically exact Eqs. J3J , © > © > 
which show that already relatively small concentration of 
wall charges dramatically reduces the self-energy trans- 
port barrier. This conclusion is in qualitative agreement 
with general statement on the role of transitional binding 
inside the channel on transit probability 10 . Our results, 
of course, can not replace powerful numerical methods 
used for description of specific biological channels 11 . 
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